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ABSTRACT 

We present spectroscopic observations of 27 active galactic nuclei (AGN) with some of the 
lowest black hole (BH) masses known. We use the high spectral resolution and small aperture 
of our Keck data, taken with the Echellette Spectrograph and Imager, to isolate the narrow-line 
regions (NLRs) of these low-mass BHs. We investigate their emission-line properties and com- 
pare them with those of AGN with higher-mass black holes. While we are unable to determine 
absolute metalliciti es, some of our obje cts plausibly represent examples of the low-metallicity 



AGN des cribed by Groves et all (120061) . based on their [N II] /Ha ratios and their consistency 
with the iKewley & Ellisonl (120081) mass-metallicity relation. We find tentative evidence for 
steeper far-UV spectral slopes in lower-mass systems. Overall, NLR emission lines in these 
low-mass AGN exhibit trends similar to those seen in AGN with higher-mass BHs, such as 
increasing blueshifts and broadening with increasing ionization potential. Additionally, we 
see evidence of an intermediate line region whose intensity correlates with L/L^dd, as seen in 
higher-mass AGN. We highlight the interesting trend that, at least in these low-mass AGN, the 
[O III] equivalent width (EW) is highest in symmetric NLR lines with no blue wing. This trend 
of increasing [O III] EW with line symmetry could be explained by a high covering factor of 
lower ionization gas in the NLR. In general, low-mass AGN preserve many well-known trends 
in the structure of the NLR, while exhibiting steeper ionizing continuum slopes and somewhat 
lower gas-phase metallicities. 
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1. Introduction 



Active galactic nuclei (AGN) with relatively small black hole (BH) masses (Mbh ^ 10 6 M ) comprise 
a demographic tha t until recently went relatively unexplore d. It is difficult to find and observe such low- 
luminosity objects ( Filippenko & Holl2003 : Barth et al ll2004 While l umino us quasars with massive central 
BHs have been studied intensively for decades (e.g. Boroson & Greenlll992h . only with the advent of large 
surveys like the Sloan Digital Sky Sur vey TSDSSilYork et allboooh has it become possible to search for large 



numbers of low-mass active galaxies ( Greene & Ho 12004 . 2007a). The number density and radiative prop- 
erties of these accreting, low-mass BHs, which we refer to as low-mass AGN, are vital to understanding the 
accretion history of the Universe. By understanding their radi ative properties, we can con strain the spectral 
energy distributions used in models of primordial seed BHs (jVolonteri & Natarajanl l2009). which may have 
prompted early galaxy formation (e.g. lBromm & Yoshidall201lL and references therei n) and were possibly 
a significant participant in the reionization of the Universe (IMilosavljevic et al.ll2009l) . In addition, inves- 
tigating these low-mass AGN is necessary to understand the overall demographics of BHs at low masses, 
including the low-mass e nd of the Mrh ~ o~* relation, which can furt her inform our understanding of galaxy 
formation and evolution (ITremaine et al.ll2002l : iGultekin et al.l l2009). 



We focus on the sample of 174 low-mass AGN selected by iGreene & Hoi (|2007al . ; hereafter GH07) 
from the SDSS. Note that these object s can be technically classifie d as narrow-line Seyfert 1 (NLS1) galax- 
ies based on their broad H/3 widths (lOsterbrock & Poggd Il985l) . These low-mass A GN and the i r hos t 



gal axies are well-stu died, with stellar velocity dispersion measurements presented in 
and 



Barth et al. 



Xiao et al 



(201 1), and ho st galaxy analysis with the Hubble Space Telescope presented by 



(2005) 



Greene et al. 



(2008) and iJiang et all (1201 lb . In addition, their spectral slopes from the UV to the X-ray h ave been used 
to assess whether the spectra l energy distribution gets harder when the BH mass is lower (IGreene & Ho 



2007b; D esroches et al. 



2009L Dong et al. in prep). In this paper we use high-resolution spectroscopy from 



the Keck telescope to specifically study the emission-line properties of 27 of these objects. The Keck spectra 
are observed with a smaller aperture than the SDSS data, using a 0.75" x 1" aperture extraction versus 3" 
fibers used in SDSS, with a ~ 22 km s -1 instrumental resolution. This improves our ability to distinguish 
emission from the narrow-line region (NLR) of the AGN compared to H II regions in the host galaxy and 
makes it possible to decompose the broad and narrow components of the permitted lines more clearly than 
is possible in the SDSS data. We make comparisons between AGN with low- and high-mass BHs, and 
investigate the metallicity of the NLR, the structure of the NLR, and the radiative properties of the low-mass 
AGN themselves. 

In Sections $2] and $3] we present our high resolution, high signal-to-noise (S/N) observations of a 
sample of 27 low-mass AGN and our fitting methodology and measurements. In Section §@J we investigate 
the gas-phase metallicities of these AGN using various emission line diagnostics in pursuit of rare, low- 
metallicity AGN. In Section §5\ we discuss the inferred far-UV continuum slope and discuss implications 
for the spectral energy distributions of accreting low-mass BHs in AGN. In Section $6j we investigate low- 
contrast, weak emission lines through use of composite spectra. In particular, we focus on high-ionization 
Fe lines and their behavior relative to more prominent, lower-ionization NLR emission lines. In Section CZJ 
we present a summary of our conclusions. 
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2. The Sample 



The galaxies studied here are a subset of those presented in iGreene & Hoi ([2004 . 2007a ). with the 
exception of J2 156+ 1103, which was drawn from the sample in lBarth et all (120051) . The low-mass candi- 
dates are AGN w ith broad Ha emission (z < 0.352) origin ally selected from the SDSS data releases one 
dYork et aljboooh and four jAdelman-McCarthy et aDbood) . So-called virial BH masses, M BH = fRv 2 /G, 
are derived indirectly using the dense gas orbiting the BH (the broad-line region or BLR) as a dynamical 
tracer. The full- width at half-maximum (FWHM ) of Ha indicates the velocity of the BLR gas and Lu a , in 
conjunction with the radius-luminosity relation (IBentz et all 120 06). gives the radius of the BLR. Fin ally, f 
is a sc aling factor that accounts for the unknown geometry of the BLR, here assumed to be / = 0.75 (INetzer 
1990). Our low-mass AGN have been selected to have an estimated Mbh < 2 x 1O 6 M . As naively ex- 
pected, t hese low-mass AGN have been found to inhabit low-luminosity (~ 1 mag below L*), late-type host 
galaxies (IGreene et al.ll2008l : LTiang et al.ll201 lh . 



The spectra in this paper are drawn from the spectra presented in iBarth et all (120051) and IXiao et al. 



(|2011r) . which are subsets of the galaxies from GH07. In those papers the focus was on the M^-a* relation 
for these systems, while here we foc us on the emissio n line properties. Note that our paper represents only 
a subset of the targets presented in IXiao et al. J201 1|) that were available when this project began. The 
general properties of each sample are the same as the parent sample in GH07. The median BH masses 
for our sample, GH07, and Xiao et al. are as follows: M BH = 3.5 x 1O 5 M , 1.3 x 1O 6 M , 9.5 x 10 5 M o , 
respectively. The median Ha luminosities are log L# Q = 40.7, 41.1, and 41.1, while the median Eddington 
ratios are L/L^^ = 0.3, 0.4, 0.3, respectively. We calculate the bolometric luminosity by applying the 
correction in GH07: L bol = 2.34 x 10 44 (L Ha /10 42 ) - 86 ergs s -1 , while L Edd = 1.26 x 10 38 (M BH /Af ). Using 
a Kolmogorov-Smirnov (KS) test, we f ind that our 27 l ow-mass AGN have statistically indistinguishable 
redshifts and luminosities compared to lXiao et all (|201ll) . (P = 0.05 and P = 0.03, respectively). Compared 
to GH07, our objects have Ha luminosities that are likely drawn from the same parent distribution (P = 0.03). 
Finally, the redshift distribution of our sample is statistically different than that of GH07 (P < 0.01), such 
that our objects are at somewhat lower redshift (median z = 0.054 compared to the median z = 0.086 for 
GH07). Note, however, that the lower redshift objects have more reliable BH mass estimates. In light of 
these overall similarities, we conclude that our objects are representative of the parent sample of low-mass 
AGN from GH07. 



Observations were taken on the Echellette Spectrograph and Imager (ESI; ISheinis et al.ll2002l) on the 
Keck telescope during observing runs in 2003, 2004, 2006, and 2008. The 0.75" slit width and 1" extraction 
of these Keck spectra potentially allow us to isolate the NLR better than in the SDSS data by excluding 
contributions from H II regions at larger radii. The observed wavelength range for the spectra is 3900 - 
10900 A over 10 echelle orders. Exposure times ranged from 900 - 1800 seconds. Note that the absolute 
flux scale of the spectra is uncertain due to slit losses and non-photometric conditions on some nights, but 
the relative flux calibration across the spectrum should be robust since the AGN and standard stars were 
observed at the parallactic angle. The average S/N of our sample at 5500 A is S/N= 26, w ith an instrumenta l 
resol ution, a ~ 22 km s -1 . The full data reduction process is described in more detail in Barth et all (120080 
and lXiao et al.1 (|201ll) . which is identical to that used here. Table[T]lists object information for our low-mass 
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AGN sample. 



3. Methodology 
3.1. Fitting Procedure 

Our investigation is focused on the properties of the nuclear gas itself, specifically the emission-line 
ratios, inferred gas kinematics, and gas-phase metallicities of the objects at the low end of the BH mass 
distribution. To this end, we want accurate measurements of the emission-line fluxes and shapes, for which 
we must remove the underlying continuum. At the low nuclear luminosities of most of our sample, starlight 
from the galaxy becomes a significant contributor to the overall continuum, and could bias the emission-line 
fluxes if not properly accounted for. In addition to the underlying galaxy, a compl ete model of the AGN 



continuum includes a nonstellar power law, typically well-fit by the form fx oc X a (Richar ds et al] |2003). 
The optical spectrum also contains emission from numerous broad Fe II transitions that form a pseudo- 
continuum underlying many other emission lines of interest (see Figure [T|). Because the components of our 
fits are coupled, our fitting process iterates between continuum and emission line fits, as described below. 

For computational simplicity, we divide the sample into two classes of objects - those where the 
AGN dominates the continuum and the galaxy is not generally measurable, and those where the galaxy 
contributes a significant fraction of the continuum (see Figures Q] and We treat the fitting of these 
AGN-dominated and galaxy-dominated spectra separately because there are strong degeneracies between 
the power-law shape, the galaxy continuum level, and extinction. We use the IRAF routine splot to mea- 
sure the equivalent width (EW) of the Ca II K A3933 A absorption feature in a user-defined window. We 
then determine an empirical division between AGN-dominated and galaxy-dominated objects based on this 
EW. After inspection of the spectra and comparison to the measured Ca II EWs, we determine that galaxy- 
dominated objects typically have Ca II K EW > 1 A, and therefore all such objects are fitted using our 
galaxy-dominated procedure^ For both AGN-dominated and galaxy-dominated fitting methods, we start 
with an initial continuum fit, th en iterate between emission lin e fitting and continuum fitting, analogous to 



the procedure to remove Fe II in lNetzer & Trakhtenbrotl (12007b . After three iterations, the Ha flux changed 



by < 0.001% for the median objects, but there were two objects where the changes were at the 10% level 
at the final iteration. Still, this 10% change is within our final errors as described below. The ratio of H/3 
fluxes between iterations also typically converged to within 0.001% of the final value, but was up to 1% for 
two objects. 



'Two objects with Ca K EW< 1A, J1434+0338 and J0931+0635, were moved to the galaxy-dominated group after inspection 
showed Mg I b absorption in their spectra. The Ca K lines were not distinctive in these cases because of lower signal-to-noise at 
the blue end of these spectra. 
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3.2. Continuum Fits 



For AGN-dominated objects, the continuum fitting windows included the following regions: 3900- 
4330, 4430-4770, 5020-5650, and 7600-7900 A. Examples of our AGN-dominated objects are shown 
in Figure Q] For these objects, we applied a continuum fit that included two power laws, Fe II pseudo- 
continuum, and reddening. There are many degeneracies in these fits between internal extincti on, Galactic 
extinc tion, flux calibration e rrors, and the intrinsic shapes of the AGN and galaxy continua. 



Barth et al. 



(|2005r) and Ixiao et al.l (120111 ) use a quadratic polynomial in addition to a power-law and galaxy continuum 
to account for the addition al terms. We have c hosen to include a single reddening law instead. We adopt the 
reddening pre scription inlCalzetti et all (2000) for star-forming galaxies, since our objects have blue, disky 
host galaxies (|Greene et all 120081) . Note that by using the Calzetti et al. reddening law we have implicitly 
assumed a dust opacity and geometry that may not be applicable to these AGN. Our original hope was 
to glean additional physically meaningful information from the fits including extinction. After comparing 
the B aimer decrements with the continuum-derived reddening, we determined that due to the degeneracies 
described above, our continuum-derived extinction terms are not robust. Therefore, we use the extinction 
term only to derive satisfactory continuum fits. 

Although our spectral region only includes 3800-8500 A, we find that two power laws are needed in 
38% of the AGN-dominated objects, where the power laws are defined as f\ oc X a . AGN often show a break 
in the slope of the continuum ar ound 5000 A, with the blue region steeper than the red, which is commonly 
fitted with a broken power law (IShang et alJl2005T) . This uptur n in the blue e nd of the continuum is not due 
to the accretion disk emission known as the little blue bump (iShieldsl Il978h . which ends around 4000 A, 
but represents a real break in the power law continuum. We cannot tell if this break is due to increasing 
galaxy continuum toward the red (IVanden Berk et alj|2001r) . because there are not enough stellar features 
to accurately model the galaxy in these AG N-dominated objects. Initial parameters for the power laws are 



based on the Ha-Lsioo relation described in lGreene & Hoi (|2005bl ). The power-law slopes were free to vary, 



but were both limited to have indices less than or equal to zero (median index a = -1.30) in order to reduce 
degeneracy with extinction. For all fits requiring a broken power law, one of the fitted power laws is flat 
with an index of zero, which dominates the red end of the spectrum, and the other power law dominates the 
blue end with an average power-law index of a = -2.26. We note again that there is still some degeneracy 
between the blue power-law slope and the level of extinction. 

We fitted t he Fe II pseudo-continuum by scaling, broadening, and shifting the empirical model from 



Borosonl (|2002l ). which is based on the prototypical strong Fe II object I Zw 1. Fitting empirical tem- 
plates from I Zw 1 is a well-established method for dealing with the Fe II in optical and UV spectra (e.g., 
Boroson & Greenlll992l : IVestergaard & Wilkesll200ll : iBorosonlbooi Leighlv & Moorel 200^ . For typical 
AGN, the I Zw 1 template fits the Fe II pseudo-conti nuum well (IKovacevic et al. While templates 



have been made using theoretical Fe II atomic models JvernerE oOO). these are limited by the the complexity 
of the Fe II ion and the lack of laboratory data on the atomic transitions. We should note that the structure 
within the Fe II pseudo-continuum diff ers from the struc ture of the I Zw 1 template in several of our galaxies, 
as seen by IKovacevic et all (120 Id!) and lHu et all (I2008bh . 



-6- 



For the galaxy-dominated spectra, we applied a continuum fit over the following wavelength regions 
to include more stellar absorption lines than the pure AGN fits: 3900-4300, 5020-5650, 5870-5920, 6150- 
6530, and 7100-7700 A. Our galaxy-dominated fits included a single power law, a model galaxy spectrum 
with an age of 10 Gyr, an d reddening. The power-law index for these objects was fi xed at a = —1.56 as i n 
Vanden Berk et al. ( 2001 ) with a starting amplitude based on the Ha-L 510 o relation ( Greene & Ho 2005bh . 
This is somewhat steeper than th e median index we find fo r our AGN-dominated objects (a = -1.30). The 
galaxy model was made using a Bruzu al & Chariot! (|2003l ) single stellar population model with an age of 
10 Gyr, where we assume a IChabrierl (12003b initial mass function, and the metallicity is assumed to be 
solar. Solar metallicity may be an over-estimate for these galaxies (see Section ® , but because of the small 
aperture of our observations, and the significant dilution from the AGN continuum, we do not have strong 
constraints on the stellar population parameters. The galaxy model is d omi nated by emission from K giants 
in the optical spectral region. The continuum fits of Barth et al. ( 2005 ) and lxiao et al] (hoill) . which utilize 
stellar templates rather than single stellar population models, also find continua dominated by K giants. The 
free parameters of this 10-Gyr galaxy model are a broadening function (a), model amplitude, and redshift . 
Note that the spectral models have an intrinsic dispersion that is larger than the galaxy spectra, and thus the 
fitted a is not physically meaningful. As above, the galaxy and power law components were multiplied by a 
reddening function, for which we varied Ay. 

For three of our galaxy-dominated objects, including J1240-0029, J2156+1 103, J0806+2419, individ- 
ual inspection revealed Fe II pseudo-continuum in the spectra, so we fitted them with the previous galaxy 
continuum method, but included the Fe II pseudo-continuum component as well. For four other galaxy- 
dominated objects (J0824+2959, J1057+4825, Jl 143+5500, and J1631+2437) an obvious A star-like com- 
ponent was necessary in the 3800-4200 A region to account for Balmer absorption. For these objects, we 
included an additional galaxy model component, much like the 10-Gyr model, b ut with an age of 1 Gyr. 
We do not find any evidence of systematic differences in host galaxy morphology ( Jiang et allboill) or line 
ratios for objects requiring a 1-Gyr galaxy model. Examples of each fitting procedure are shown in Figure |2] 



3.3. Emission Line Fits 

Once the continuum fitting is complete, we have line-only spectra to fit. The emission line-fitting 
procedure was the same across all continuum types. We fit narrow and broad components to the Balmer 
lines, from Ha to He, narrow and broad He I lines A5876, A6678, and A7065, the [S II] AA6717, 6730 
doublet, the [N II] AA6548, 6584 doublet, [O I] lines A6300 and A6364, [O III] lines AA5007, 4959 and 
A4363, narrow and broad He II A4686, and N III A4640. The fitted emission lines are illustrated in Figure [3] 
and are summarized in Table |2j where we list the line species, rest wavelengths, and the components used in 
the fit. We present derived quantities of our low-mass AGNs in Tabled and the emission-line measurements 
for all of our objects in Tables [4] and [5] Note that the fluxes presented have only been corrected for Galactic 
reddening. Widths of multi-component lines are calculated by constructing the total fit of the line and 
measuring the FWHM of this reconstruction. 
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For each object, we fit the NLR emission lines using sets of Gauss-Hermite functions, much like 



Salviander et all ([2007). The expression for a set of Gauss-Hermite functions is as follows: 



F(x) = Ae- xl '\\ +h 3 Mx) + h 4 f 4 (x)] (1) 

Mx) = -^(2V2x 3 - 3V2x) (2) 
v6 

/ 4 W = ^=(4x 4 -12x 2 + 3) (3) 
V24 



where 

x=(X-X ceat )/a. (4) 

When using a single set of Gauss-Hermite functions, each emission line is fit with five parameters: line 
centroid, width (a), flux, h 3 , and h 4 . The h 3 parameter quantifies the asymmetry of the line and h 4 measures 
the kurtosis, or boxiness, of the line. Therefore, we can in principle quantitatively compare NLR line shapes 
among objects. This works well for NLR emission with the exception of [O III], where 74% of our objects 
are not well-described by a single set of Gauss-Hermite functions, as described below. In general, we base 
the NLR line shapes on a fit to the [S II] doublet because it is in a relatively clean region o f the spectrum, is 



not blended with o ther lines, and is usually strong enough to constrain the line shape (e.g. iHo et al.lll997bl : 



Greene & Hd 12004) . Unless otherwise noted, all NLR lines are fixed to have the same velocity width, h^, 
and h 4 as the [S II] fit. The wavelengths for each species are fixed to their relative laboratory ratios. The 
flux ratio of the [N II] doublet is fixed to its laboratory value of 2.96. We also fix the [O III] A4363 line to 
have the same shape and number of components as the [O III] A5007 fit, such that its flux is the only free 
parameter. In two cases, J0806+2419 and 1534+0408, the narrow H/3 line was so weak that we imposed a 
limit that H/3 could not exceed Ho:nlr/3. 

For the [O III] A5007 line, we try an initial fit with one Gauss-Hermite function for each line, with 
starting values taken from the [S II] fit, but all five parameters free. As described above, [O III] lines are 
often complex in AGN and can have strong wavelength shifts and/or asymmetries, so we allow the code to 
fit the lines with two summed Gauss-Hermite functions, and again allow the code to determine if there is a 
statistically significant improvement with the additional component. We find that 74% of our spectra require 
two components for the [O III] fit, and of those, 85% had a blueshifted second component. 

While we tried fitting the BLR lines with Gauss-Hermite functions as well, we foun d that more than 
one G auss-Hermite component is required to model the broad lines. This is in contrast with lsalviander et aD 



d2007D . where they were able to model Ha and H/3 with single sets of Gauss-Hermite functions. Our 
observations have much higher resolution and S/N than the SDSS spectra they used, and thus we have 
robust measurements of two component BLRs. Therefore, we chose to fit the broad lines with sums of 
Gaussians, because the interpretation of summed Gauss-Hermite functions becomes no more intuitive than 
sums of Gaussians, and in fact, has more degeneracies among the fit parameters. 
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For the broad Ha, we allow the code the freedom to use up to four Gaussians an d to pick the mos t 
effective fit for the number of free parameters included. Using a criterion inspired by IHao et all ([2005), 
the code determines the benefits of adding each additional Gaussian to the fit by comparing the number of 
additional model parameters to the improvement in x 2 - We apply this same procedure to the broad H/3, with 
up to three possible summed Gaussians. Examples of our Ha and H/3 fits are shown in Figure |4] We then 
subtract the fit of the Ha complex from the line spectrum, in order to fit the He I lines. 



3.4. Error Analysis 

To estimate uncertainties on our emission line fits, we generate Monte Carlo simulations of each galaxy 
spectrum. To construct the model, we take the measured values from the final continuum and emission line 
fits and reconstruct the spectrum, to which we add Gaussian random noise using the measured error array. 
For each object, we construct 100 such fake spectra and fit them using our continuum and emission line 
fitting prescription. The final uncertainty is defined by the distribution of measurements of the artificial 
spectra. Specifically, we use the width encapsulating 68% of the results. 

Of course, we also incur systematic errors in our continuum fitting. To estimate their magnitude, we 
explore the consequences of removing components from the most complicated fits: those requiring a young 
(1-Gyr old) galaxy component. We explore fits without the 1-Gyr continuum component and with or without 
an Fell pseudocontinuum component. Even in this rather extreme case, the flux ratio of the narrow Ha/H/3 
lines typically only varies by 10% or less across all the different fits, although one object showed a factor of 
two increase with the addition of the 1-Gyr galaxy component. 

We also compa re our fitted measurements of the narrow and broad components of Ha to those made 
bv lxiaoetal]j201lh on many of the same objects, and find that 93% of our measured widths are consistent 



to within la and 75% have fluxes consistent to within 20% of the measured value. This provides a valuable 
additional measure of systematic errors and is consistent with our measured uncertainties. 



4. Gas-Phase Metallicity 



Because AGN are some of the most distant observable object s, metallicities in AGN ar e used as a 
possible tracer of the history of chemical evolution in the universe (IHamann & Ferlandlll999l) . While in- 
vestigation of high-redshift AGN typically probes the BLR metallicity, here we focus entirely on the metal- 
licity of the NLR, which should be more representative of the gas-phase metallicity of the central region 
of the gal axy. AGN typically have super-solar metallicities, as measured with diagnostic emission lines in 
the NLR. broves et all <|2006h found that out of a sample of ~ 23,000 AGN from SDSS, only 40 showed 



sub-solar metallicities. To some degree, this result is biased by the difficulty of finding low-mass BHs in 
late-type, star-forming galaxies with sub-solar gas-phase met allicities. Only a handful of low-metailicity , 



dwarf galaxies which harbor AGN have been detected to date (IKraemer et al.l l 1999: 



Thus, we still do not know how intrinsically rare these low-mass BHs are (e.g. 



Izotov & Thuan 



Greene & Ho 



2008). 



2007a). We 
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suspect that when AGN are found in low-mass galaxies, the gas-phase metallicitie s will be low and in accord 
with the mass-metallicity relation of inactive galaxies (e.g. iTremonti et all 120041) . We test that supposition 
directly here using the NLR emission lines. 

We investigate the gas-phase metallicities of the objects in our low-mass, low-luminosity sample using 
several rest-frame optical NLR emission-line diagnostics. Because nitrogen is a secondary element, as the 
abundances of metals relative to hydrogen increase, the abundance of nitrogen relative to other metals also 
increases. As a result, the ratio of nitrogen to hydrogen is a good metallici ty diagnostic, and in particular. 



rNni/i 

1986b. 



[Nil] /Ha has traditionally been used in research on AGN and other nebulae ( Pagel & Edmunds! 1981 ; Evans 
The loc us of low-metallicity objects is easily seen in the well-known Baldwin, Phillips, Terlevich 



(BPT) diagram (Baldwin et al.lll98ll) . For a given [O IIIJ/H/3, AGN that have a lower [N II] /Ha are likely 
at lower relative metallicity. 

For the sake of comparison to typical AGN, we compare our sample of low-mass AGN with a sub- 
set of the SDSS DR7 galaxy sample at z ~ 0.1. Spectral measurements for the SDSS subset have been 
taken from the catalog compiled by the Max-Planck-Institute for Astrophysics/Johns Hopkins University 



Brinchmann et al. 



(MPA/JHU) groujffi. The fitting process used to construct the catalog is described in 
(120041) and ITremonti et al] (120041) . Note that the catalog selects narrow-line AGN. We select only objects 
whose Ha, H/3, [O III], and [N II] emission lines had errors < 10% of the flux of the lines, which results in 
a total of 147,816 objects, plotted in grayscale in Figured Star-forming galaxies lie to the left of the dotted 
line at log [N II] /Ha < -0.2, and AGN are in the plume extending to the upper right. Two of our galaxies, 
Jl 143+5500 and J2156+1 103, fall into the star-forming region as denoted by the Kauffman et al. dashed 
line, indicating contamination from H II regions within the host galaxies. While the presence of broad emis- 
sion lines in our objects confirms the presence of an AGN in these two galaxies, they must have substantial 
contributions from star formation within the host galaxy. Compared to local AGN in DR7, our sample of 
low-mass AGN extend to lower [N II] /Ha than the SDSS AGN, which seems to in dicate that the low-mass 
AGN include objects at relatively lower metallicities than the bulk of SDSS AGN. [H useman n et al. (1201 ll) 



also see this trend with the large sample of low-mass SDSS AGN presented in iGreene & Hoi (12007 



M). 



We want to determine absolute metallicities for our low-mass AGN to confirm whether they do rep- 
resent the appa rently rare, sub-sol ar A GN. However , it is di fficult to define metallicities for AGN on an 
absolute scale. iDopita et al.1 (120021') and iGroves et all d2004ah have developed photoionization models that 
include the photoionization structure and emission from dusty, radiation pressure-dominated gas. To this 
end, they employ MA PPINGS III, a photoionizatio n and shock code whose development is detailed in 
Dopita et all (119820 and lsutherland & Dopital J 1993b . MAPPINGS III considers patches of gas to have the 
whole range of ionizati on structure , there by producing emission which contributes to virtually all of the lines 



observed in an NLR. IDopita et al.l (|2002) argue that by causing the radiation pressure to dominate the pres- 
sure gradient of the gas, the ir models can predic t emission line ratios that are consistent with the observed 
values seen in real systems dGroves et al.ll2004bh . Note that MAPPINGS III reports a total metallicity, not a 
gas-phase metallicity, so the MAPPINGS III metallicities may be systematically higher than other literature 



'http://www.mpa-garching.mpg.de/SDSS/DR7/ 
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values. For a more extensive overview of photoionization modeling methods, see also lGrovesi (120071) . 



We estimate metallicities in our objects by comparing the [O III] /H/3 vs. [N II] /Ha measurements for 
our sample of AGN to the predictions from the MAPPINGS III grid of models. For our sample, we are able 
to estimate the density in the NLR using the [S II] doublet (100 <n e < 1300, see Table and to constrain 
the shape of the AGN continuum based on previous observations of a for AGN (a = -1.56; Vanden Berk 
et al. 2001). This allows us to utilize models that are appropriate for our particular sample. These dusty, 
radiation pres sure-dominated mode ls are publicly available within the library of grids provided in the IDL 
code ITERA (IGroves & All en 2010). In Figure [6l we show our low-mass AGN measurements compared to 
four grids of models, which represent the full variety of physical parameters we see in our objects, including 
n e = 100 cm" 3 and 1000 cm" 3 , and a = -1.4 and -1.7. While there are some regions of the grids that can be 
described by several models, especially at high log U and steep indices (a = -1.7), generally the metallicity 
increases with increasing [N II] /Ha. The metallicity estimates derived from the MAPPINGS III code can 
vary by up to a factor of three in metallicity, depending on the range of a considered. MAPPINGS III 
metallicity estimates for individual objects can range from sub-solar to several times solar depending on the 
assumed spectral shape. Our metallicity estimates are presented in Table [3] There are also several objects 
that lie at the upper edge of the models, which may reflect limitations of the models. 

To do a relative comparison, we investigate trends between our low-mass AGN and the SDSS sample of 
AGN using the MAPPINGS III models. To derive the metallicity for each of our galaxies, we consider where 
an object lies on a particular MAPPINGS III grid, and interpolate between the two nearest metallicity tracks 
to determine a metallicity (Z) for that object (see Figure©. For each object, we estimate the metallicity from 
grids with the nearest density, n e = 100 cm -3 or n e = 1000 cm -3 . The average change in metallicity between 
models with n e = 100 cm" 3 and n e = 1000 cm" 3 at fixed a = -1.4 is only 0.02 Z , so density differences 
are negligible. However, the difference in metallicity for a = -1.4 versus a = -1.7 is 0.3 Z Q , and since the 
errors in our slopes are large, we incorporate this uncertainty into our results. We estimate the metallicity 
using both a = -1.4 and a = -1.7, but exclude estimates for the a = -1.7, Z = 4 Z & model, because it is 
degenerate with the other models in the region of interest (see Figure 0. This gives a range of Z for each 
object presented as black bars on Figure [7] 

To highlight the difficulties in deriving absolute metallicities, we consider the particular case of NGC 
4395, a well-studied, low-mass AGN. We ha ve an ESI spectrum o f NGC 4395, for which extensive measure- 
ments and Cloudy modeling are published (IKraemer et al.ll 19991) . and we find that our observed emission- 
line ratios agree with the previous measurements [log ([N II]/Ha) ~ -0.5, log ([O IIIJ/H/3) ~ 1.0]. 
Kraemer et al.1 (119991) estimate a metallicity of Z ~ 0.25 Z Q , which is consistent with estimates made from 
the H II region s. Since then, the d efinition of solar metallicity has been revised because of new oxygen 
measurements (|Asplund et al.ll2004l ). such that solar metallicity is now defined to be a fa ctor of two lower 
than previously thought. Therefore the estimate of metallicity from IKraemer et al. Nl999l) actually supports 
a metallicity for NGC 4395 of 0.5 Z . To estimate a metallicity from the MAPPINGS III grids, we need 
to determine the density and power-law index for NGC 4395. Using the [S II ] A67 17, 6731 double t, we 



estimate that the density in the NLR of NGC 4395 is about 1300 cm" 3 , and IKraemer et all dl999) also 



constrains a„ 



-1.7 for NGC 4395. Using these values, we estimate the metallicity of NGC 4395 to be 
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Z ~ 1.5 Z from the MAPPINGS III models that have n= 1000 cm" 3 and a = -1.7 (Figure[6d]>. There is a 
significant zero-point offset between the two sets of photoionization models, as has been well described in 
other contexts. Again, we note that the absolute scale of the gas-phase metallicities remain highly uncer- 
tain. For instance, one reason the MAPPINGS III values appear high is that the depleted metals are counted 
(i.e., we calculate the total, rather than gas-phase metallicity). Other factors, such as the photoionization 
modeling methodology, are likely at play as well. In the fo llowing, we use me tallicity estimates from the 
MAPPINGS III models so that we can compare directly with lGroves et all (|200q) . but focus only on relative 
trends between ob j ects r ather than absolute metallicities. Our approach is analogous to that advocated by 
Kewley & Ellisonl (2008) for star-forming galaxies. 



For comparison to AGN with m ore massive BHs, we selected AGN from the SDSS DR7 galaxy sample 
described above using the cuts from Iho et all (|l997ah for Seyferts ([S II] /Ha > 0.4 and [O IIIJ/H/3 > 3). 
This results in 3747 galaxies, which are binned according to the stellar mass reported by the MPA/JHU 
group in increments of 0.5 dex between 9.0 < log M*/M & < 12.0. For each stellar mass bin, we estimate 
Z/Zq from the MAPPINGS III grids using the median [N H]/Ha and [O II1J/H/3 value for that bin (red bars 
in Figure |7] denote the range of metallicities that result from the range of a as above). We also estimate the 
range of Z for the central lu of the objects in each bin in the same way (grey bars in Figure [7]), except for the 
highest- and lowest-mass bins which had fewer than 40 objects. Because these SDSS AGN cover the same 
range of densities as our low-mass AGN (300-1500 cm -3 ), we estimate their Z from the same MAPPINGS 
III grids as described above. 

We present our relative metallicity comparison in the context of the M*-Z relation in Figure [7] The 
M*-Z relation is only defined for star-forming galaxies and has never been measured in AGN host galaxies 
before to o ur knowledge . For s everal of our low-mass AGN, we derive a stellar mass from the Hubble Space 
Telescope (Green e etalj 120081) /-band luminosities and a mass-to-light ratio based on the SDSS g-r color 
following Bell et al. ( 2003 ). These stellar masses and our metallicity estimates are presented in Table [3^1 To 
expand the low-mass end of the sample, we inclu de NGC 4395 for compari son, with metallicities estimated 
from the MAPPINGS III grids. We also plot the Kewley & Ellisonl (120081) M*-Z relation for star-forming 
galaxies that is based on [O III] and [N II] (referred to as PP04 03N2), which is the closest analog to our 
method but was derived for a star-forming spectral energy distribution. For those low-mass AGN with 
10.0 < log M*/M & < 10.5, which includes most of our objects, we find that the mean (Z/Zq) = 0.11, which 
is somewhat higher than the Z predicted by the M*-Z relation for star-forming galaxies shown by the solid 
curve. In comparison, the metallicities of our objects agree with the SDSS AGNs at the same mass. 

Even without modeling, our low-mass AGN include objects that have a lower ratio of [N II] /Ha relative 
to the SDSS AGN (Figure |5]), which suggests lower metallicities. This confirms our suspicion that since our 
host galaxies extend to lower mass than the SDSS galaxies, the low-mass AGN include lower-Z objects. 
There appears to be a systematic difference between the active and inactive galaxies, in the sense that the 
active ones appear to have higher average metallicities at fixed stellar mass (Figure [7]). Either we are seeing 



3 We also calculate stellar masses using a fixed mass-to-light ratio of 0.6 in the /-band (motivated by Greene et al, 20081) to 
estimate the lower limit of the stellar masses. With such a low mass-to-light ratio the masses would decrease by 0.4 dex. 
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systematic differences due to methodology (i.e. different ionizing spectra) or it is a real effect. Since the 
NLR is compact c ompared to the interstellar medium of the galaxy as a whole at these luminosities (e.g. 
Schmitt et al.1l2003h . we are probing the hi gh-metallicity central r egion of the galaxy. Given well-known 

Ih 



metallicity gradients in spiral galaxies (e.g. lHenry & W orthev 1999), we may be slightly overestimating the 
overall gas-phase metallicity in the active galaxies. 

Our investigation shows that determining absolute metallicities for AGN is a very complex problem, 
and as suggested in other contexts, we can only trust relative metallicity information derived using the 
same methodology across samples. While the AGN in our sample have similar NLR metallicities as the 
SDSS AGN at fixed stellar mass, our low-mass AGN include objects at somewhat lower metallicities than 
SDSS AGN, at least in part because they occupy lower-mas s host galaxies. We estimate metallicities using 



[N II] /Ha and photoionization models. 



Barth et al. 



present the Seyfert 2 counterparts to our low- 



mass Type 1 AGN, which have overlapping [N II] /Ha and [O IIIJ/H/3 values and therefore similarly likely 
include low-Z objects as well. Because of the uncertainties from th e models we cannot say for sure whether 
our objects are the rare, sub-solar metallicity objects addressed by iGroves et al.1 (|2006T) . Conversely, when 
directly compared to the metallicity predictions from the M*-Z relation for star-forming galaxies, AGN are 
consistent with obeying a comparable mass-metallicity relation as inactive galaxies, although the scatter is 
large. 



4.1. Metallicity Systematics 

As we will see in §5\ there is tentative evidence for a steeper far-UV continuum slope in these systems 
of a = -2.0. If we assume a steeper slope of a = -2.0 for the photoionization models, this would only 
minimally impact the metallicity estimate of objects with estimates Z < 1 Z Q , causing Z to increase by 
~ 0.2. A steeper slope of a = -2.0 could potentially increase the measured metallicities by up to a factor 
of two for objects with Z > 1 Z©, but this does not change our previous c onclusions th a t som e of these 



low-mass AGN potentially represent the rare, low-metallicity AGN sought in lGroves et al.l ((2006). 



Line emission from star formation in the host galaxy could also potentially affect our metallicity es- 
timates. To estimate the magnitude of the changes, we take two representative AGN with metallicities of 
solar and twice solar. We then add varying l evels of star formation, assuming gas at the same metallicity 
and deriving the appropriate line ratios from Pettini & Pagel ( 20041) . For instance, an AGN with a solar 



gas phase metallicity will have [O IHJ/H/3 « 6, while the star forming regions will have [O III]/H/3 « 0.2. 
In the extreme case that 75% of the B aimer emission comes from star formation, the observed line ratio 
would be [O III]/H/3 = 1.7. For solar metallicity, we underestimate the metallicity estimates only by 10-30% 
for a 25-75% contribution from star formation to the Balmer line flux. At twice solar, 50-75% contribu- 
tion from star formation leads to a factor of 2-3 underestimate in metallicity. However, only two sources, 
Jl 143+5500 and J2156+1103, lie far enough from the AGN locus to be consistent with this level of con- 
tamination (Fig. [5]). Thus, we conclude that contributions from star formation do not change our metallicity 
distribution significantly. 
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5. Spectral Energy Distribution 



In theory, the lower masses of these BHs should correspond to a hotter accretion disk temperature and 
thus measurable differences in the spectral energy distributions (SEDs) of our sample from those of more 
typical AGN. Testing this expectation will help to bridge the gap in BH behavior and accretion theory from 
stellar mass BHs to the largest SMBHs. In addition, the nature of the accretion disk surrounding these low- 
mass AGN can inform our understanding of the primordial inter mediate-mass BHs that fo rmed within the 
first galaxies and participated in the reionization of the Universe (|Milosavljevic et alJl2009f ). We know that 



the SE Ds of accretion disks around stellar mass BHs are harder and flatter than those of quasars (Frank et al. 



2002). In fa ct, the peak blackb ody temperature for disks in AGN with low-mass BHs may also move into the 
X-ra v band (IDone et alJl2012h . We have see n hints of SEP changes in this sample i n previous work, both in 
a ox JPesroches et al.ll2009r) and in EW M icroom et alJ l20o3: breene & Holl2005bl) ; While most of our ob 



jects have been stud ied at radio ( Greene et al.ll2006h and X-ray JGreene & Hol2007bl : Desroches et al. 1 120091 : 
Miniutti et al.ll2009l) wavelengths, we have no direct constraints on the far-UV continuum. We can gain an 
indirect handle on the accretion disk shape, however, since those photons are responsible for photoionizing 
the emission line gas that we observe in the optical. 

Specifically, we can infer the slope of the continuum in t he far UV by contr asting the He II and narrow 
H/3 components. As described in IPenston & Fosburyl (|1978h . since both H/3 and He II are recombination 
lines, the ratio of their relative strengths should yield a measure of the relative intensity of the continuum at 
912 A and 228 A, the ionization edge for He II. Assuming the ionizing source is well-described by a power 
law of the form /„ oc u avv , the intensities should behave as follows: 



^4686 
Ih/3 



= 1.99x4 auv 



(5) 



As the BH mass decreases at fixed Eddington ratio, we would expect to move closer to the peak of the 
big blue bump, and thus to observe shallower auv as a function of mass. Of course, m any other factors ma y 
be at play, including differing Eddington ratio distributions or different reddening (e.g. lBonning et al. 2007). 
The deduced UV power-law slopes for our objects range from -4 < auv < -0.4, with a media n value of 
-1.9. To compare to a sample of typical AGN, we consider the low-redshift AGN presented by iHo & Kim 
(|2009h . which have Mbh ~ 10 6 - 10 8 M Q and a median absolute magnitude of Mb ~ -19 derived from their 
H/3 luminosities. Because their observations were taken on the Magellan 6.5 m Clay Telescope, they have 
data of similar quality to ours. They measured many emission lines, including He II and H/3. We use their 
measured line ratios to calculate auv for their 94 objects, and present the distribution of the two samples in 
Figure [8] The comparison sample has a range from -2.5 < auv < 0.3, with a median slope of -1.4. A KS 
test of the two distributions results in P < 0.001, meaning that there is less than 0.1% probability that they 
are drawn from the same parent sample. We find that the distributions are statistically different whether we 
include the three low-mass AGN with upper limits in our distribution or not. 

Given this difference in distributions, and the significant shift in the median, we also look for a cor- 
relation between Mbh and auv- We find a mild correlation, with a correlation coefficient r = 0.18 for the 
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range of BH masses included here, from 5.5 < log Mbh/M < 10. A least-squa r es fit to the data yields the 
following relationship: log Mbh °c 0.10 auv- Previous work (iShang et al.ll2005l : Davis et alj|2007r) investi- 
gate possible correlations between the UV continuum slope and Mbh and do not find compelling evidence 
for a trend between UV slope and Mbh- We probe a bluer region of the UV continuum and extend to lower 
M B h than those authors, and we see a tantalizing hint for a correlation between M B h and auv- On the other 
hand, the sense of the trend is counter to our naive expectations for lower-mass BHs. It would be useful to 
compile a sample including our low-mass B Hs with matched Eddin gton ratios and examine both the X-ray 
and UV spectral slopes at the same time (e.g. iDesroches et al.l l2009). 



6. Composite Spectra 



An effective way to highlight general spectral characteristics is to make composite spectra. Other au 
thors i n AGN research have use d composite spectra to investigate trends in AGN across a wide range of red 
shifts (IVanden Berk et al. boo lb . to probe weak emission lines in high luminosity AGN (Francis et al. 1991) 
and w eak emission and absorption lines in z = 2 - 3 star-forming galaxies that house AGN ( Hainline et al 
201 ID . One advantage of a composite spectrum with our high spectral resolution is that it gives us the 
opportunity to investigate low-contrast features that would otherwise be difficult to detect, including faint 
emission lines, line wings, and intermediate-line regions. Some of our objects exhibit high-ionization, for- 
bidden Fe lines, such as [Fe VII] A6087, [Fe X] A6374, and [Fe XI] A7892. These coronal lines are thought 
to originate in the inner edge of the NLR, where densities are low enough t o allow forbidden line emission, 
but the ionization parameter is high enough to enable ions like Fe +9 to exist. iMullaney et al.l (|2009) compare 
the line shapes and shifts of the centroid of these lines to [O III]. They find that the high-ionization lines are 
at similar offsets, and therefore similar velocities, to the blue wing of \Q III], which may originate in out- 
flowing material (|De Robertis & Osterbrocklll984l : Icecil et al.ll2002l : belbord et al.ll2009h . With the luxury 
of very high spectral resolution and S/N, we can both investigate the trends highlighted in previous works in 
more detail and determine whether they hold for these low-mass, low-luminosity systems. 



6.1. Composite of Entire Sample 



We construct our composite spectra by interpolating all the spectra to the same rest-wavelength grid 
using [S II], normalizing them to the rest-frame AGN continuum at 5600 A, and then taking the median flux 
value of contributi ng objects at each wavele ngth. This method is similar to the median composite construc- 
tion presented in IVanden Berk et all (|2001h . To highlight the weak features, we construct our composites 
from the continuum-subtracted spectra. In Figure |9l we present the composite spectrum of the original spec- 
tra including all 27 low-mass AGN, which highlights the many weak NLR emission that we observe, and in 
Table 5 we present our measurements of the many NLR emission lines in this composite. 

In our continuum-subtracted composite, which highlights emission line features, we fi rst investigate 
broad Ha and H/3, to compare the BLRs in our low-mass AGN to those of more massive AGN. Greene & Ho 
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(12005m investigate the relationship between the FWHM of Ha and H/3 for a sample of over 200 AGN from 
the SDSS, which were selected to have high S/N and low galaxy contamination. They find that H/3 is 
commonly broader than Ha, and fit a relation between the two. Within our composite spectrum, we also 
find that the wings of the H/3 line are broader than those of Ha. This trend can be seen in our individual fits 
as well, with a median FWHM of Ha of 840 km s" 1 and a median FWHM of H/3 of 1195 km s" 1 . In the 
accepted model of the BLR, where the emission lines are primarily broade ned by their Keplerian veloc ities, 
H/3, being broader, is emitted from gas interior to the region emitting Ha (|Osterbrock & Ferlandll2006l . and 
references therein). 

Turning to the NLR, we investigate the line shapes of some of the strongest NLR e mission lines. 



includ ing [O III], [S II] A6731, and the high-ionization [Fe VII], to see if the results from iMullaney et al. 



d2009h hold for AGN of lower mass (see Figure [lOab . While [S II] and [O III] have very similar profiles on 
the red si de, [O III] shows subst antial excess flux on the blue side. It turns out that [Fe VII], much like the 
profiles in Mul laney et al. lies mostly under the blue wing of [O III]. The common velocity structure 

of the blue wing and the high ionization lines supports the idea that the high-ionization Fe lines are emitted 
in an outflow from the inner face of the dusty torus, and appears to hold true even for our smaller scale 
low-mass AGN. 



We also present multiple high-ionization Fe lines in Figure I10b[ which all have ionization poten- 
tial > 100 eV. In the composite of all 27 low-mass AGN, the [Fe VII] line is blueshifted by ~ 75 km s -1 
and the [Fe X] line blueshifted by ~ 125 km s -1 relative to [S II]. The [Fe X] line is also broader than the 
[Fe VII] line. Unfortunately, the composite shows only a marginal detection in the [Fe XI] line, but the 
[Fe VII] and [Fe X] lines suggest that as the ionization level of these Fe lines increases, they get broader 
and m ore blueshifted. This pattern has been seen in high- ionization, forbidden lines in AGN by ma ny au- 
thors (brandilll978l booke et al.lll97d : IPenston et al]ll984h . In comparison to belbord etall (T20091) . who 
investigate 63 AGN from SDSS which are selected to have high-ionization forbidden lines, the increasing 
FWHM and blueshift of high-ionization Fe lines in our low-mass AGN most closely resemble their NSL1 
and Seyfert 1.5 subsets. Increasing FWHM and blueshift at high ionization potential have also been iden- 
tified using near-infrared, high - ioniza tion, forbidden lines of various elements, including S, Si, Fe, Ca, and 
Al, by iRodnguez-Ardila et al.1 (|201ll ). One would expect this behavior if, as the high-ionization, forbidden 
lines increase in ionization potential, they ar e emitted from closer to the nucleus than oth er NLR lines, and 
are nearer the launching point of an outflow (|Heckman et al.ll 198 ll : [Ward & Morrisll 19841 ). A possible alter- 
native explanation is that the high-ionization emission results from inflowing material on the far side of the 
central engine, which is still blueshifted from our point of view; however, the higher ionization lines would 
still be spatially closer to the nucleus with their stronger blueshifts. 

Because the high-ionization Fe lines suggest a possible relationship between ionization potential and 
line widths and shifts, w e investigate several lower-ioniz ation, narrow lines to look for trends with ionization 
potential, as discussed in lOsterbrock & Mathews! (1198 6) and references therein. We fit the composite spectra 
with our emission-line fitting procedure as described above. In addition, we measure line widths and shifts of 
weaker lines presented here using a single Gaussian fit to each line in the continuum-subtracted composite 
spectrum, because there is not enough signal in the high-ionization lines to make a more complicated fit 
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informative. We do this analysis on the composite spectrum rather than individual spectra because the high- 
ionization lines are weak enough that they are difficult to measure in individual spectra and benefit from the 
combined signal in the composite. In our individual spectra, only about 20% of our objects have detectable 
[Fe X] or [Fe XI], despite 16 of 27 showing [Fe VII]. 

In Figures lllal and II lbl we show the measured line widths and shifts for a selection of NLR emis- 
sion lines of varying ionization states in the composite spectrum. In accordance with the trend we see in 
Figure UObl as the ionization potential of the lines increases, the lines become broader (correlation coeffi- 
cient of r = 0.88) and more blueshifted (r = -0.56). Again, this is consistent with an ionization structure 
of the NLR in which the lines with higher ionization potential are emitted closer to the central engine and 
thus are broader. We postulate that the blueshift indicates an outflowing component on these small scales 
( Komossa et all2008 : Muller-Sanchez et"alll2011 ). Note that nuclear outflows are not typically aligned with 
the host galaxy in any systematic way. In larger samples, no correlation is see n between the presence of 
a blue wing and the inclination of the host galaxy (e.g. Greene & Ho l2005ah . and radio jets, which in- 
dicat e the direction of n uclear outflows, seem to be randomly oriented with respect to their host galaxies 
(e.g. IKinney et al. hooch . Another commonly discussed correlation is that between F WHM and critical den- 
sity. W e also find a weak trend between critical density and FWHM, as discussed by IFilippenko & Halpern 
(119841) and others. However, with a correlation coefficient of r = 0.54, this correlation is weaker than that 
with ionization potential, at least in this sample. 



6.2. Division by Physical Properties 

While the previous composites are useful for evaluating characteristics of our entire sample, we can 
also investigate how our spectra change in relation to specific physical or spectral characteristics. To this 
end, we divide our sample into several subsets. We investigate many characteristics in this manner, including 
luminosity, Eddington ratio, Mbh, FWHM of Ha, presence or absence of a blue wing in [O III], shift of the 
centroid of [O III], I13 and h\ of the NLR fits, and NLR density determined from the [S II] line ratio. In 
most cases no interesting trends were seen, or they were redundant with those that we show here based on 
Eddington ratio, luminosity, and the presence of a blue wing in [O III]. 

For each characteristic, we divide the entire sample of 27 objects into two subsets as described be- 
low and construct a continuum-subtracted composite spectrum of that subset . For th e luminosity sub- 



sets, we derive luminosities from Lh« using the formalism in iGreene & Hd (1200 5b) to calculate Lb i- 
We then divide the sample in half at the median luminosity (log Lt, i = 40.5 erg s -1 , with a full range of 
39.3 erg s -1 < log Lt, i < 41.5 erg s -1 ). To investigate Eddington ratio, we bifurcate the sample about 
the median Eddington ratio (log L/Lgdd = -0.4, with a full range of -1.8 < log L/Lsdd < 0.3). To evalu- 
ate spectra with or without a blue wing in [O III], we visually inspect the individual continuum-subtracted 
spectra to determine the presence or absence of a blue wing in [O III] in each object and then form a com- 
posite of objects with a blue wing and another composite for those without. Dividing by visual inspection is 
necessary because of the difficulty of parameterizing blue asymmetry with measurements from the two com- 
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ponent [O III] fits. For each of these six composite spectra, we investigate emission line behavior among the 
NLR emission lines shown in Figure [TTal We present the velocity dispersion in the NLR (<r) for each subset 
in Table |7J where we subtracted the instrumental resolution (<7i nst = 22 km s" 1 ) in quadrature. These values 
are averages over all of the NLR emission lines shown in Figure [TTal or averages of the high-ionization Fe 
lines, where the ionization potential > 100 eV. 

Investigation of these subsets yields interesting results for both the BLR and the NLR. We find that 
the BLR emission lines differ significantly between the high-L/LEdd and the low-L/Lndd composites. In 
Figure[l2l we show the Ha and H/3 regions for both composites. While the wings of the lines (v > 1200 km 
s -1 ) are similar, the high-L/LEdd composite shows excess flux at intermediate velocities in both the broad 
Ha and H/3 emission lines. The difference in intermediate BLR component is unique to the L/L^dd subset. 
The intermediate- width component in the high-L/LEdd composite has a FWHM= 840 km s -1 , compared 
to the width of the broad component at FWHM= 2,400 km s -1 . The intermediate line component is still 
broader than the average NLR width, shown in Table 6, and also broader than the high-io nization Fe lines. 



Our discovery of an intermediate line region supports the conclusion of iHu et al.1 (|2008ah that changes in 
Eddington ratio among objects can account for differences in an intermediate component of the Balmer lines. 
Increased intermediate-line emi ssion in high-L/LRdd AGN could bia s their derived BH masses, causing 
underestimates of the true M B h ( ICollin et alJbood : 1 Wang & Weilbo09l) . 



We also find that NLR emission-line properties are dependent on luminosity, L/LEdd, and the presence 
of a blue wing in [O III]. When investigating the effects of luminosity, we find that the high-Lboi composite 
has broader NLR emission lines than both the low-Lhni composite and the genera l comp osite of our whole 
sample. This is in keeping with past findings Jphillips et all 1 19831 : lwhittlelll985L fl992l) . where higher lu- 
minosities result in broader line wi dths, commonl y thought to be due to the inferred larger gravitational 
potential of the bulge of the galaxy (IHo et al.ll2003h . The high-L/LEdd subset also possesses similarly broad 
NLR emission lines compared to the low-L/LEdd subset. Given the small dynamic range in Mbh represented 
here, it is not surprising that we see similar trends in Lt, i and L/L^dd- 

There are well-known correlations between the line widths and line shifts of [O III] as a function of 
Eddington ratio , such that high-L/LEdd objects tend to have broader [O III] with a more blueshifted peak 
(|Borosonl 120051). This is often interpreted in the context of disk winds or outflows that are more prevalent 
at high L/L E dd dDe Robertis & Osterbrocklll984[) . Blueshifted [O III] clouds at high velocities (up to 3200 
km s -1 ) are seen in NGC 1 068, and are associated with outflow s induced by the jet and/or radiation pressure 
from the AGN itself (e.g. ICecil et al.ll2002l : iGreene et al.ll2012h . Such radial motion might be expected to 
reduce the covering factor of NLR gas and lower the EW of lines like [O III]. The well-known Eigenvector 
1 (EV1) ties these trends together (IBoroson & Greenlll992r) . At the high L/L^dd end of EV1, the ratio of 
[O III] /H/3 drops and the incidence of blue asymmetries in H/3 incr eases. The [0 III] EW may drop in high 
EV1 objects as well, but the data remain inconclusive on this point (|Boroso n 2005; Ludwig et al.ll2009r ). In 
a related tren d, <tnt.r/<7* may increase with increasing L/L^dd due to increasing nonvirial motions caused 



<7* seems to disappear altogether (IGreene et al. 2009). 



by outflows (IGreene & Hdl2005al:lHoll2009r). At the h ighest luminosities, the correlation between o"nlr and 
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Examining our low-mass AGN alone, we do not find a consistent story connecting the NLR to Edding- 
ton ratio. Turnin g first to the ratio of ctnlr/o-*, there does not appear to be any correlation with L/LEdd 
(IXiao et all 120111) . Unfortunately, the s tellar velocity d ispersions are often overestimated due to contam- 
ination from the galaxy disk (see also LTiang et al.ll201 11) . which complicates our interpretation. Secondly, 
we investigate the presence of a blue wing in [O III] and find no correlation with L/L E dd (r = 0.06). If, as 
we suggest above, blue wings are indicative of an outflowing co mponent, then we do not see evidence for 
stronger outflows at higher Eddington ratios, as we might expect (|Proga et al . 2000). Furthermore, if L/L^dd 
were tied to outflows, we might expect to find lower EW narrow-lines at higher L/L^dd- Instead, we see 
higher EW NLR lines in the higher L/L^dd composite. Perhaps our dynamic range in L/L E dd (or M B h or 
luminosity) is too narrow to truly evaluate the interdependence of outflows and accretion rate. 

On the other hand, we do see evidence that blue wings are associated with outflows. We might expect 
that in objects wit h strong outflows, the gas would be more disturbed, leading to a lower covering factor and 
thus a lower EW (baskin & Laorll2005l : lLudwig et al.ll2009t) . We do find that the EW of [O III] in our objects 
is lower, relative to the AGN power-law continuum, when there is a blue wing in [O III] (Figure [T3l) . The 
EWs of other low-ionization NLR emission lines are also lower (Figure [T4l). Furthermore, if we divide the 
sample into those with and without blue asymmetry in the low-ionization lines (e.g., the I13 measurement of 
[S II]) then we also see low EW narrow-lines in the blue-asymmetric subsample. Thus, we do see indirect 
evidence that blue wings indicate outflowing material (as we suggested above based on the velocity of the 
high ionization lines, 36.11) . Interestingly enough, all the low-ionization narrow lines, not just [O III], appear 
to behave similarly, both in terms of EWs and asymmetric lines. 

From our sample alone, we see evidence for varying levels of disturbance in the NLR, but we see 
no direct tie with L/L E dd or luminosity. Given the limited dynamic range in luminosity, and the large 
uncertainties on Eddington ratio, it is hard to dra w strong conclusions from these data alone. Perhaps the 
absolute luminosities of the sample are too low (IGreene & H o 2005b) or the uncertainties on our L/L^dd 
estimates are too large. 



7. Conclusions 



We present observations of a sample of 27 low-mass AGN (10 4 M Q < M BH < 2 x 10 6 M ), observed 
with the Echellette Spectr ograph and Imager on the Keck Telescope. Large samples of low-mass AGN have 
not existed until recently ( Greene & Ho 20041) . and so we compare their emission-line properties to those of 
well-studied higher-mass AGN. 

We investigate the NLR metallicities of these objects using emission line ratios, particularly [N II] /Ha 
and [O DXJ/H/3. Our sample includes objects with weaker [N II] /Ha than higher-mass AGN for a given 
[O III]/H/3, which implies that those objects have lower metallicities. Thus we see so me galaxies with sim- 
ilar metallicities to the rare AGN with sub-solar metall icity sought by Groves et al.1 (120061) . yet we cannot 



determine their metallicities on an absolute scale (e.g. iKewley & Ellison! I2008f) 7 Additionally, we do see 



weak evidence for a correlation between galaxy stellar mass and gas-phase metallicity in these systems. 
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Most likely, the lack of low-metallicity AGN high lighted bvlGroves et al.1 (120061) is attributable to the diffi- 
culties of finding AGN in low-mass galaxies (e.g. iGreene & Holl2007ah rather than extra enrichment of the 
NLR by the AGN itself. 

We examine the continuum properties of the accretion disk for our low-mass AGN. By using the emis- 
sion from the recombination lines [He II] and H/3, we infer the slope of the far-UV continuum from 200- 
900A, near the peak of the big blue bump. We find that the far-UV slopes in our low-mass AGN are steeper 
than those of low-r edshift, higher-mass AGN. While we know that ou r low-mass AGN are very radio-quiet 
(IGreene et al.l l2006) and have flat X-ray slopes (IDesroches et al. I l2009h compared to typical AGN, they have 
quite similar optical properties and unexpectedly, somewhat steeper UV slopes. These steeper inferred UV 
continua imply that changi ng continuum shape could explain the inverse Baldwin effect seen between H/3 
FWHM and luminosity (|Croom et al.ll2002l : breene & Holl2005bh . 



Despite tentative evidence for different SEDs, the NLR structure is also similar between low-mass and 
typical AGN. Using co mposite spectra, we ar e able to measure the widths and velocities of high-ionization 
Fe lines. As found by iMullaney et al. J2009h . the [Fe VII] line exhibits a similar width and blue shift as 
the blue wing of [O III], pointing to a common physical origin for both transitions in a radially flowing 
component. As seen in previous work, we find that the width of NLR lines and their blueshift correlates 
with the ionization potential of the line. We see no evidence for a dramatic change in the NLR structure in 
this mass and luminosity regime. 

Making composite spectra in subsets of luminosity, Eddington ratio, and presence or absence of a blue 
wing in [O III], we find first that the high-luminosity composite has much broader NLR emission lines than 
the low-luminosity counterpart within our sample, as has been seen in more-massive AGN. We also find that 
the subset showing the presence of a blue wing in [O III] exhibits weak NLR emission. We posit that in the 
objects with a blue wing, the ou tflow drives gas away from the central region, leading to a lower covering 
fraction (e.g. iLudwig et al.ll2009r) . However, contrary to our expectations from typical AGN, the presence 
of a blue wing does not correlate with high Eddington ratio, possibly due to our small dynamic range in 
this sample. Lastly, we find that the high Eddington ratio composite has excess emission at int ermediate 



veloci ties in the Balmer lines. This corresponds to the ILR identified by iHu et al.1 (|2008ah and Izhu et al. 



(|2009f ). This ILR emission could potentially bias Mbh estimates for AGN at high L /L^m by decreasin g the 
measured FWH M of the Balmer lines, thus underestimating Mbh, as pointed out by lCollin et all (I2006h and 
Zhu et all J2009h . 



Overall, the low-mass AGN studied here seem to behave in much the same way as more massive AGN, 
as traced by optical emission lines. They are likely to have relatively low gas-phase metallicities in the NLR 
than their larger brethren, and may have steeper far-UV continuum slopes, but the structure and organization 
of their emission line regions seems to be largely similar. 
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Table 1 . Low-Mass AGN Observations 



Object 


Flag 


z 


S/N 


Fit Class 


SDSS J010712. 03+140844. 9 


GH01 


A AT7A 

0.0770 


17 


AGN 


SDSS J024912. 86-081525. 6 


GH02 


0.0296 


18 


AGN 


SDSS J0325 15.59+003408.4 


GH03 


0.1021 


12 


AGN 


SDSS J080629. 80+241955. 6 




0.0415 


31 


GAL 


SDSS J080907. 58+441641.4 




0.0540 


24 


GAL 


C C C T A O 1 CCA T> . OCA/' /I A A 

SDSS J081550.23+250640.9 




0.0726 


20 


GAL 


SDSS J082443. 28+295923. 5 




0.0254 


54 


GAL 


C 1 C 1 O TA O O A 1 O . CAA/'CO O 

SDSS J082912. 67+500652. 3 


GH04 


0.0435 


37 


GAL 


r~i r-v n ri TAA ,1 O 1 A 1 O . s f\ A C C C\ 1 

SDSS J094310. 12+604559.1 


GH05 


0.0745 


25 


AGN 


C C 1 C T1A1 1 A O A A , A AO A AO T 

SDSS J 101 108.40+002908.7 


GH06 


0.1001 


12 


AGN 


SDSS J105755.66+482502.0 




0.0732 


26 


GAL 


SDSS J110501.97+594103.6 




0.0336 


57 


GAL 


SDSS Jl 14343.76+550019.3 




0.0271 


26 


GAL 


SDSS J122342.81+581446.1 




0.0146 


34 


GAL 


SDSS J124035.81-002919.4 


GH10 


0.0810 


31 


GAL 


SDSS J125055.28-015556.6 


GH11 


0.0815 


20 


GAL 


SDSS J13 13 10. 12+051942.1 




0.0488 


24 


GAL 


SDSS J141234.67-003500.0 


GH13 


0.1271 


17 


AGN 


SDSS J143450.62+033842.5 


GH14 


0.0283 


20 


GAL 


SDSS J153425. 59+040806.7 




0.0395 


9 


AGN 


SDSS J162636.40+350242.0 




0.0341 


38 


GAL 


SDSS J163159.59+243740.2 




0.0436 


43 


GAL 


SDSS J170246.09+602818.9 


GH16 


0.0691 


20 


GAL 


SDSS J172759. 15+542147.0 


GH17 


0.0995 


14 


AGN 


SDSS J215658.30+1 10343.1 




0.1080 


39 


GAL 


SDSS J232159.06+000738.8 


GH18 


0.1838 


12 


GAL 


SDSS J233837. 10-0028 10.3 


GH19 


0.0356 


29 


GAL 



Note. — Low Mas s AGN Sampl e Obse rvations. Col. (1): SDSS 
Designation. Col (2): Greene &Hol i 2004 ) object designations. Col 
(3): Redshift. Col. (4): Median S/N per pixel at 5500 A. Col (5): Fit 
classification, as described in Section §[3] 




Fig. 1. — Upper panel. Emission line-subtracted spectrum of J0107+1408. Emission lines are shown with 
dotted lines, the total continuum fit is shown in red, including the combined reddened power laws as the 
dashed blue line, and the reddened Fe II component is shown in purple. Middle panel. Continuum fit of 
J1412-0035 with the same convention as above. Lower panel. Continuum fit of J0943+6045 with same 
conventions as above. 
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Fig. 2. — Upper panel. Emission line-subtracted spectrum of Jl 105+5941. Emission lines are shown with 
dotted lines, the total continuum fit shown in red, including the reddened power law as the dashed blue line, 
and the reddened 10-Gyr galaxy component in green. Middle panel. Spectrum of J0806+2419 with the 
same convention as above, where the additional reddened Fe II component is shown in purple. Lower panel. 
Spectrum of J0824+2959 with same conventions as above, where the additional reddened 1-Gyr galaxy 
model is shown in yellow. 
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Fig. 3. — Upper panel. Original spectrum of J0829+5006 with total continuum fit shown in red, including a 
single power law as the dashed-blue line, and 10-Gyr galaxy continuum in green, which have been multiplied 
by a reddening law. Lower panel. Zoom in of continuum-subtracted spectrum to enhance the view of weaker 
emission lines. All fitted emission lines are labeled. 
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Table 2. Fitted Emission Lines 



Emission Line 


Wavelength (A) 


Possible components 


He I, narrow 


7065, 6678, 5876 


1 G-H Poly each 


He I, broad 


7065, 6678, 5876 


1 Gaussian each 


[S II] 


6730, 6717 


1 G-H Poly each 


[Nil] 


6584, 6548 


1 G-H Poly each 


Ha - He, narrow 


6563, 4861, 4341, 4102, 3970 


1 G-H Poly each 


Ha, broad 


6563 


< 4 Gaussians 


[0 I] 


6364, 6300 


1 G-H Poly each 


[OIII] 


5007, 4959, 4363 


< 2 G-H Polys each 


H/3, broad 


4861 


< 3 Gaussians 


He II, narrow 


4686 


1 G-H Poly 


He II, broad 


4686 


< 2 Gaussians 


NIII 


4640 


1 Gaussian 


H7 - He, broad 


4341, 4102, 3970 


1 Gaussian each 



Note. — Col. (1): Emission Line Species. Col (2): Rest wavelength (A). 
Col (3): Possible components in fit, where G-H Poly is a set of Gauss- 
Hermite functions. 
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Table 3. Derived Measurements 



Object 


Lot 


Mbh/Mq Log L/L Edd 


Log Af ga i/M 


n e (cm 3 ) 


Z/Zq 


JO 107+ 1408 




6.03 


-0.40 




360 


2.0±0.1 


J0249-0815 




5.32 


-0.66 




1010 


2.0±0.3 


J0325+0034 




5.98 


-0.57 


10.05 


230 


1.5±0.5 


J0806+2419 




^ CK 
J.yJ 


— U. / 1 


10.24 


100 


1.6±0.2 


J0809+4416 




5.94 


-0.61 


10.04 


140 


1.8±0.4 


J0815+2506 




5.55 


-0.50 


10.24 


170 


1 .8 ± 0.5 


J0824+2959 




5.33 


-0.93 


10.16 


590 


1.5±0.2 


J0829+5006 




5.76 


-0.12 




270 


2.3 ±0.7 


J0943+6045 




5.76 


0.02 




190 


1.9±0.2 


J101 1+0029 




6.18 


-0.56 




210 


0.8 ±0.2 


J1057+4825 




5.68 


-0.30 


10.12 


60 


0.9±0.1 


Jl 105+5941 




6.89 


-1.81 


10.68 


440 


1.5 ±0.3 


Jl 143+5500 




5.97 


-1.25 




500 


0.5±0.1 


J1223+5814 




6.04 


-0.93 


9.65 


850 


1.6±0.2 


J1240-0029 




5.93 


-0.09 




400 


1.6 ±0.4 


J1250-0155 




6.80 


-1.25 




160 


0.9 ±0.2 


J1313+0519 




5.32 


-0.40 


10.17 


340 


2.9±1.0 


J1412-0035 




6.15 


0.11 




1030 


1.8±0.2 


J1434+0338 




5.65 


-0.68 




240 


2.5 ±0.6 


J1534+0408 




4.79 


-0.21 


9.36 


360 


2.4 ±0.6 


J1626+3501 




5.54 


-0.39 


10.32 


190 


2.6±1.4 


J1631+2437 




5.33 


0.02 


9.92 


160 


2.2 ±0.8 


J1702+6028 




6.16 


-0.31 




10 


0.9±0.1 


J1727+5421 




5.68 


-0.14 




290 


1.6±0.6 


J2156+1103 




6.45 


-0.65 




230 


0.6±0.1 


J2321+0007 




6.56 


-0.78 




290 


2.2 ±0.8 


J2338-0028 




5.92 


-1.14 




540 


2.0±0.1 


Note. — Col. (1): Object name. Col. (2): Log M bh /Mq, from Xiao et al. (2011). Col. 
(3): Log Ljcz/Z-Edd calculated from our Hciblr measurements. Col. (4): Log Mt/Mp, for 
obiects with reported /-band luminosities fromlGreene et alj J2008) (see ifil. Col. (5V. Density, 
calculated from the [S II] doublet. Col. (6): Z/Zq, as determined in 



Table 4. Emission Line Measurements: Ha Region 



V ) UJ CV^ L 




[S II] A6717 




rc TTi \A7^1 

|_J -L-LJ y\\J 13 1 
















J 1 J( n P) 


FWHM 


h 3 




f /ffFf/31 

J / JK n P> 


J I JK n P) 


FWHM 
i vv nivi 


Av 


J 1 J\ n P> 


J/J( n P> 


FWHM 

1 VV JTLLVl 


Av 


J0107+1408 


0.63 ±0.04 


90 ±20 


-0.03 


-0.03 


0.57 ±0.04 


2.6 ±0.3 


120 ±40 





0.63 ±0.04 


50.6 ±2.6 


960 ±400 


-10 


J0249-0815 


0.66 ±0.04 


110±40 


0.01 


0.05 


0.75 ±0.05 


3.8±0.3 


130 ±30 


23 


0.66 ±0.04 


18.3 ± 1.0 


790 ± 250 


-68 


J0325+0034 


0.36 ±0.03 


120 ±30 


-0.08 


0.01 


0.30 ±0.03 


2.9 ±0.3 


150 ±20 





0.36 ±0.03 


21.5±0.9 


1070 ±190 


-229 


J0806+2419 


0.36 ±0.35 


70 ± 120 


-0.13 


0.10 


0.29 ±0.28 


6.0±5.8 


220 ± 140 


-2 


0.68 ±0.65 


104.7 ± 100.3 


1120±210 


31 


J0809+4416 


0.67 ±0.02 


160 ±10 


-0.01 


-0.04 


0.53 ±0.02 


3.4±0.1 


180 ±20 


12 


0.67 ±0.02 


9.7 ±0.2 


1020 ±80 


81 


J0815+2506 


0.56 ±0.03 


130 ±40 


-0.10 


0.04 


0.45 ±0.03 


2.9 ±0.2 


150 ±40 





0.56 ±0.03 


15.9±0.7 


840 ± 1 10 


-22 


J0824+2959 


0.64 ±0.02 


270 ±10 


0.05 


-0.13 


0.64 ±0.02 


6.1 ±0.2 


310± 10 


24 


0.64 ±0.02 


3.3±0.1 


1240 ±20 


-236 


J0829+5006 


0.48 ±0.01 


150±10 


0.02 


-0.05 


0.41 ±0.01 


2.7±0.1 


160 ± 10 





0.48 ±0.01 


13.5±0.2 


690 ±30 


-33 


J0943+6045 


>0.29±0.10 


140 ±30 


-0.08 


-0.02 


> 0.23 ±0.08 


> 1.4±0.5 


150 ±40 


59 


>0.45±0.16 


> 48.1 ±16.1 


660 ± 420 


-22 


J101 1+0029 


0.20 ±0.02 


100 ±20 


-0.01 


0.04 


0.16 ±0.02 


2.7 ±0.2 


210 ±20 


81 


0.20 ±0.02 


28.8 ±0.8 


1220 ±250 


46 


J1057+4825 


0.60 ±0.02 


80 ±10 


-0.04 


0.06 


0.44 ±0.02 


3.8±0.1 


70 ±10 


23 


0.60 ±0.02 


3.6±0.3 


640 ± 140 


-45 


Jl 105+5941 


0.80 ±0.02 


300 ±20 


0.07 


0.06 


0.74 ±0.02 


4.1 ±0.1 


350 ±20 





0.80 ±0.02 


3.4±0.1 


4240 ±780 


-125 


Jl 143+5500 


0.61 ±0.01 


70 ±10 


-0.01 


0.04 


0.43 ±0.01 


3.3±0.1 


70 ± 10 


23 


0.21 ±0.01 


1.1 ±0.1 


1290 ±70 


-171, 


J1223+5814 


0.51 ±0.01 


100±10 


-0.04 


0.00 


0.55 ±0.01 


4.6±0.1 


140 ±20 


12 


0.51 ±0.01 


12.2±0.2 


1150± 150 


-22>J 


J 1240-0029 


0.51 ±0.01 


1 AC\ -4- 1 n 


U.LJl 


0.03 


0.47 ±0.01 


4.2±0.1 


160 ± 10 





0.51 ±0.01 


10.1 ±0.1 


840 ±80 


-183" 
1 


J1250-0155 


0.71 ±0.01 


110±10 


-0.14 


0.07 


0.56 ±0.01 


3.5 ±0.1 


110± 10 


-2 


0.36±0.01 


1.9 ±0.1 


2260 ±160 


11 


J1313+0519 


0.80 ±0.02 


150±10 


-0.01 


0.03 


0.71 ±0.02 


3.5 ±0.1 


180 ±10 





0.80 ±0.02 


6.6 ±0.2 


720 ± 180 


-56 


J1412-0035 


>0.30±0.11 


150±30 


-0.15 


0.04 


> 0.35 ±0.12 


> 2.1±0.8 


180 ±40 


-215 


>0.53±0.18 


> 57.5 ±19.2 


670 ±370 


-23 


J1434+0338 


0.64 ±0.03 


120 ±30 


-0.05 


0.04 


0.53 ±0.03 


3.4 ±0.2 


120 ±20 


23 


0.71 ±0.02 


14.4 ±0.4 


860 ± 190 


-68 


J1534+0408 


0.55 ±0.07 


70 ±20 


0.01 


0.00 


0.49 ± 0.07 


3.0 ±0.4 


100 ± 10 


10 


0.55 ±0.07 


34.4 ±3.5 


480 ± 130 


24 


J1626+3501 


0.40 ±0.02 


70 ±20 


-0.02 


-0.03 


0.32 ±0.01 


3.6±0.1 


90 ± 10 


23 


0.68 ±0.02 


24.2 ±0.5 


750 ±40 


-45 


J1631+2437 


0.56 ±0.02 


180 ±20 


0.02 


0.07 


0.45 ±0.02 


3.4±0.1 


210 ±20 





0.56 ±0.02 


10.5 ±0.2 


600 ± 140 


-332 


J1702+6028 


0.47 ± 0.02 


180±30 


-0.06 


-0.15 


0.33 ±0.02 


4.5 ±0.2 


260 ±20 


-22 


0.48 ±0.02 


31.5±0.9 


1040 ±160 


-321 


J 1727+5421 


0.41 ±0.02 


150 ±60 


-0.08 


-0.01 


0.35 ±0.01 


3.1±0.2 


200 ±70 





0.41 ±0.02 


12.3±0.3 


810 ±690 


12 


J2156+1103 


0.18±0.01 


160 ±30 


-0.08 


-0.01 


0.15±0.01 


2.8±0.1 


190 ±20 


11 


0.24 ±0.01 


6.7±0.1 


1540 ±230 


120 


J2321+0007 


0.65 ±0.06 


180±110 


-0.05 


-0.01 


0.56 ±0.07 


3.6 ±0.3 


200 ±80 


12 


0.65 ±0.05 


27.5 ±1.4 


1670 ±1020 


-114 


J2338-0028 


1.22±0.01 


120±10 


0.01 


0.07 


1.19±0.01 


3.8±0.1 


130 ± 10 


12 


1.22±0.01 


2.6±0.1 


1490 ±80 


-33 



Note. — Col. (1): Object name. Cols. (2-6): Fitting parameters for [S II] A6717. Col. (7): Flux of [S II] A6731, where all other parameters are fixed to [S II] A6717. Cols. (8-10): Fitting parameters 
for Hcinlr, where and A4 are fixed to [S II] A6717, and Av is the centroid shift compared to [S II]. Col. (11): Flux of [N II], where all other parameters are fixed. Cols. (12-14): Parametric 
measurements for total Hqblr fit. FWHMs and Av are in km s , and fluxes are relative to H/3nlr, which can be found in Table[5] Fluxes have been corrected for Galactic extinction using values from 



ISchlegel et alj jl998l) and the reddening law of lO'Donnelll jl994l) . Upper and lower limits are three sigma estimates. 



to 



Table 5. Emission Line Measurements: H/3 Region 



\J UJ CC I 


[0 III] A5007 










He II A4686JVLR 










///(H/3) 


FWHM 


Av 


t (rip) 


FWFTM 
r vv nivi 


Av 


///(H/3) 


FWHM 


Av 


f 1 fiWR\ 
J 1 J\ n P) 


FWHM 
r vv nivi 


Av 


J0107+1408 


12.9 ±0.7 


150± 10 


7 


0.14±0.01 


160 ±30 


30 


0.29 ±0.05 


150 ±30 


1.2 


14.3 ±0.8 


1380 ±330 


18 


J0249-0815 


11.1 ±0.6 


110± 10 


30 


0.15±0.01 


160 ±30 


41 


0.33 ±0.05 


150 ±30 


1.3 


4.7 ±0.3 


1190 ±380 


-15 


J0325+0034 


6.3 ±0.3 


100 ± 10 


7 


0.13±0.01 


400 ±20 


-A 


0.16±0.04 


380 ±20 


12.7 


3.8±0.2 


3120 ±260 


375 


J0806+2419 


12.0 ± 11.5 


180 ±20 


-24 


0.08 ±0.06 


230 ±1250 


24 


0.13±0.29 


220 ±1250 


11.5 


38.4±36.9 


1530 ±670 


14 


J0809+4416 


3.2±0.1 


160 ± 10 


-4 


0.56 ±0.01 


190 ±20 


41 


0.10±0.02 


180 ±20 


1.2 


2.8±0.1 


1550 ±280 


179 


J08 15+2506 


3.6 ±0.2 


130 ±20 


18 


0.11 ±0.01 


180 ±40 


30 


0.14±0.03 


170 ±40 


1.3 


4.4 ±0.2 


1410 ±440 


30 


J0824+2959 


16.8±0.3 


360 ± 10 


38 


2.83 ±0.03 


340 ± 10 


32 


0.15±0.01 


320 ± 10 


-0.2 


0.8±0.1 


730 ±90 


-32 


J0829+5006 


6.3 ±0.1 


180 ± 10 


18 


1.56±0.01 


200 ± 10 


18 


0.17±0.01 


190 ± 10 


1.3 


3.0±0.1 


1040 ±40 


18 


J0943+6045 


>4.4±1.5 


190 ±20 


-19 


< 0.27 ±0.08 


160 ±40 


47 


>0.53±0.18 


600 ±40 


26.9 


> 14.6±4.9 


760 ± 240 


47 


J101 1+0029 


4.4 ±0.2 


160 ± 10 


30 


0.19 ±0.01 


360 ±30 


110 


0.43 ±0.03 


340 ±30 


1.3 


7.3 ±0.3 


1690 ±290 


87 


J1057+4825 


1.7 ±0.1 


300 ±50 


-96 


0.34 ±0.01 


80 ±20 


41 


0.01 ±0.01 


80 ±20 


12.9 


0.6±1.5 


720 ±1090 


-4 


Jl 105+5941 


2.5±0.1 


400 ± 30 


-4 


6. 10 ±0.03 


370 ±20 


30 


0.03 ±0.01 


350 ±20 


1.2 


0.8±0.1 


6740 ±1980 


594 


Jl 143+5500 


1.5 ±0.1 


70 ±10 


30 


1.07 ±0.01 


80 ±10 


53 


0.01 ±0.01 


70 ±10 


1.4 


0.1 ±0.1 


460 ± 10 


317 


J1223+5814 


12.3 ±0.2 


100 ± 10 


18 


1.40±0.01 


160 ± 10 


30 


0.31 ±0.01 


150± 10 


1.3 


1.9±0.1 


1730 ±400 


41 


J1240-0029 


/ it U. 1 


1 4- in 


18 


1.99 ±0.01 


170 ±10 


30 


U. 14 ± U.UI 


1 An -i- 1 n 


1.2 


2.4 ±0.1 


1130 ±100 


-130 


J1250-0155 


4.0±0.1 


150± 10 


5 


1.79 ±0.01 


120 ± 10 


26 


0.07 ±0.01 


110± 10 


4.8 


0.6±0.1 


10170 ±6670 


563 


J1313+0519 


5.0±0.1 


160 ± 10 


-A 


0.48 ±0.01 


200 ± 10 


30 


0.14±0.01 


190 ± 10 


1.3 


1.3 ±0.1 


11 10 ±260 


-4 


J1412-0035 


>5.4±1.9 


190 ± 10 


-21 


< 0.16 ±0.04 


120 ±30 


-288 


> 0.07 ±0.05 


110 ±30 


4.8 


> 19.4±6.5 


950 ±400 


8 


J1434+0338 


6.1±0.2 


140 ±10 


-4 


0.31 ±0.01 


150 ±20 


41 


0.16±0.03 


140 ±20 


1.2 


2.8±0.1 


1220 ±270 


64 


J1534+0408 


6.7 ±0.7 


80 ±10 


18 


0.05 ±0.01 


200 ±20 


40 


0.13±0.08 


190 ± 10 


8.6 


9.5 ±1.0 


500 ±210 


49 


J 1626+3501 


3.7±0.1 


90 ±130 


18 


0.36 ±0.01 


100 ± 10 


53 


0.12±0.02 


100 ± 10 


1.3 


3.9 ±0.2 


910± 150 


7 


J1631+2437 


6.2±0.1 


210± 10 


18 


2.13±0.03 


220 ± 20 


30 


0.13±0.02 


210 ±20 


1.2 


3.2±0.1 


400 ±50 


-326 


J1702+6028 


4.4 ±0.2 


200 ±10 


-27 


0.76 ±0.02 


420 ±30 


-39 


0.13±0.03 


400 ±30 


24.2 


5.6 ±0.2 


3890 ±940 


-15 


J1727+5421 


5.3 ±0.2 


170 ±20 


18 


0.64 ±0.01 


250 ±40 


30 


0.24 ±0.02 


240 ±40 


1.4 


2.8±0.1 


1300 ±690 


53 


J2156+1103 


0.5±0.1 


440 ± 170 


-201 


1.82±0.01 


220 ± 20 


45 


< 0.31 ±0.11 


50 ±20 




1.0±0.2 


4260 ±470 


^08 


J2321+0007 


4.5 ±0.3 


210 ±40 


-15 


0.08 ±0.01 


250 ±70 


30 


0.21 ±0.04 


230 ± 70 


1.2 


8.3 ±0.5 


3930 ±1200 


76 


J2338-0028 


13.6±0.1 


120 ± 10 


30 


1.29 ±0.01 


140 ±10 


41 


0.37 ±0.01 


130 ± 10 


1.2 


0.5±0.1 


1200 ±470 


7 



Note. — Col. (1): Object name. Cols. (2-4): Parametric measurements for total [O III] A5007 fit. Cols. (5-7): Fitting parameters for H/3/vlr- Note that F(H/3) is in units of 
10 -15 ergs s~' cm~ 2 A -1 , and hj and /14 are fixed to [S II] A6717. Cols. (8-10): Fitting parameters for [He II]jvlr, hj and /14 are fixed to [S II] A6717. Cols. (11-13): Parametric mea- 
surements for total H^blr fit. FWHM s and Av are in km s - ', and fluxes are relative to H/3jvLfi. Fluxes have been corrected for Galactic extinction using values from lSchlegel et al] jl998h and 
the reddening law of O'Donnell 1 1994). Upper and lower limits are three sigma estimates. 
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4600 4800 5000 6500 6600 6700 

Rest Wavelength (A) 

Fig. 4. — Left column. H/3 fits, where dotted purple lines show narrow components of H/3, He II and [O III]. 
Blue, cyan, and green lines are broad components of H/3, He II, and N III. Right column. Ha fits, where 
dotted purple lines are [N II] and Honlr fits, and blue, cyan, and green show components of Hoblr- In the 
insets we zoom in on the [S II] fit. For both columns, the panels show (from top to bottom) J0824+2959, 
J0806+2419, Jl 105+5941, and J0829+5006. 
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Fig. 5. — BPT diagram of SDSS DR7 galaxy sample in grayscale contours, and our low-mass AGN as 
crosses, with low-mass AGN NGC 4395 and POX 52 for comparison. Also plotted is the theoretical 



Kauff mann et al. (2003) divider (dashed) between starburst galaxies (below) and AGN (above), and the 



dotted 



Ho et al. 



(1 1997a!) lines, which separate Seyferts (upper right) from LINERs (lower right). 
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Fig. 6. — MAPPINGS III models used to determine the metallicities in our sample, where Z = 0.25,0.5, 1,2, 
and 4 Z Q , from left to right, and log U ranges from -4 to from bottom to top (Top row: n= 100 cm -3 , 
Bottom row: n= 1000 cm -3 , Left column: a = -1.4, Right column: a = -1.7). Using these grids, our sample 
shows a range of metallicity estimates, with most objects having super-solar metallicities. 
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Fig. 7. — Host galaxy stellar mass vs. NLR metallicities for our sample (Fig. © as solid black bars. The 
range in Z represents the effect of a reasonable range in assumed ionizing slope in the modeling. Solid red 
bars are obscured AGN from the SDSS DR7 with metallicities deri ved in the same way. The la variation 
over each mass bin is shown in grey. The solid curve represents the iKewley & Ellisonl ([2008) stellar mass- 
metallicity relation based on [O III] and [N II], which is bracketed by dotted lines showing the most extreme 
relations they present (KK04 is uppermost and P05 is lowermost). 
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uv 



Fig. 8. — Distribution of derived UV slope, from He II to H/3 ratio, for the GH objects (solid) and the Ho & 
Kim sample (dotted). 
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Fig. 9. — Continuum-subtracted composite spectrum of low-mass AGN showing weak emission line features 
listed in Table [6] 
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Table 6. Composite NLR Emission Lines 



Line A FWHM f/f H/3 Ion. Pot. 



[Ne III] 


3868.7 


165 


0.42 


63.5 


[FeV] 


3889.0 


100 


0.07 


75.0 


[Ne III] 


3967.2 


153 


0.03 


63.5 


He 


3970.2 


112 


0.06 


13.6 


[S II] 


4068.9 


150 


0.04 


23.3 


H8 


4101.9 


112 


0.16 


13.6 


H7 


4340.6 


112 


0.26 


13.6 


[OIII] 


4362.9 


158 


0.13 


54.9 


He II 


4685.9 


199 


0.19 


54.4 


H/3 


4861.1 


208 


1.00 


13.6 


[OIII] 


5006.7 


157 


4.35 


54.9 


[FeVII] 


5720.3 


272 


0.03 


125.0 


He I 


5875.8 


115 


0.09 


24.5 


[FeVII] 


6085.6 


223 


0.05 


125.0 


[OI] 


6300.3 


112 


0.19 


13.6 


[OI] 


6364.1 


112 


0.07 


13.6 


[FeX] 


6371.4 


459 


0.05 


262.1 


[Nil] 


6548.1 


162 


0.38 


29.6 


Ha 


6563.0 


111 


2.02 


13.6 


He I 


6678.2 


115 


0.02 


24.5 


[S II] 


6716.7 


112 


0.48 


23.3 


[S II] 


6731.1 


112 


0.44 


23.3 


He I 


7065.3 


115 


0.03 


24.5 


[Arlll] 


7135.8 


105 


0.09 


40.7 


[Arlll] 


7751.6 


106 


0.01 


40.7 


[Oil] 


7320.1 


160 


0.04 


35.1 


[Oil] 


7330.1 


126 


0.02 


35.1 


[FeXI] 


7890.5 


332 


0.05 


290.2 


[OI] 


8445.5 


255 


0.06 


13.6 


[S III] 


9069.1 


104 


0.32 


34.8 


[S III] 


9531.0 


125 


0.63 


34.8 



Note. — Narrow emission-line measurements of 
continuum-subtracted composite spectrum of low-mass 
AGN. Col. (1): Line ID. Col. (2): Measured centroid in 
A. Col. (3): FWHM in km s" 1 . Col. (4): Flux relative 
to narrow H/3. Col. (5): Ionization potential in eV. 
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Fig. 10.— left panel: Line profiles of [O III] in black (solid), [S II] A6731 in red (dashed), and [Fe VII] 
A6087 in blue (dot-dashed) in velocity space. The line fluxes have been scaled for easier visibility, right 
panel: Line profiles of [Fe VII] A6087 in black (solid), [Fe X] A6375 in red (dot-dashed), and [Fe XI] A7892 
in blue (dashed) in velocity space. In both panels, v = is with respect to [S II]. 
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Fig. 11. — Top: Line widths (FWHM) plotted for NLR forbidden lines from the continuum-subtracted 
composite of low-mass AGN. Bottom: Line shifts (relative to the centroid predicted from [S II]) for narrow 
lines from the continuum-subtracted composite of low-mass AGN. Negative shifts are blueward. 
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Table 7. Subset Properties 







Ave. 


a (km s ') 


Composite 


No. of objects 


All NLR 


High-ion Fe 


Whole sample 


27 


70 


135 


High L bo i 


14 


98 


188 


Low L bo i 


13 


53 


98 


High L bo i/L E dd 


14 


90 


175 


Low Lboi/ L E dd 


13 


55 


102 


Includes [O III] Wing 


17 


85 


136 


No [O III] Wing 


10 


66 


131 



Note. — NLR emission-line properties for subset composites. Av- 
erages are taken over all measured NLR emission lines, or over high- 
ionization Fe lines (ionization potential > 100 eV) as noted. 
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Fig. 12. — Composites of the Ha {left panel) and H/3 regions {right panel) for high L/Ledd (solid) and low 
L/LEdd (dotted) objects. Notice the excess ILR in the high L/L^dd composite for both lines. 
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Fig. 13. — Composite spectra with a blue wing in [0 III] (solid) and without (dotted). In order to compare 
EWs, we plot each continuum-subtracted composite plus the average power-law AGN continuum of objects 
in that composite, where the power law is normalized to 5600 A. Note that symmetric [O III] lines have 
much higher EW. 
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Fig. 14. — Logarithmic NLR emission line EW ratios vs. ionization potentials. We compare the EWs in the 
composites with and without a blue wing. 



